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Abstract 

We do parametric calculations to elucidate multi-terminal electron transport properties through a molec- 
ular system where a single phenalenyl molecule is attached to scmi-inhnite one-dimensional metallic 
leads. A formalism based on the Green's function technique is used for the calculations while the model 
is described by tight-binding Hamiltonian. We explore the transport properties in terms of conductance, 
reflection probability as well as current-voltage characteristic. The most significant feature we articulate 
is that all these characteristics are very sensitive to the locations where the leads are connected and also 
the molecule-to-lead coupling strengths. The presence of other leads also has a remarkable effect on these 
transport properties. We study these phenomena for two-, three- and four-terminal molecular systems. 
Our numerical study may be utilized in designing tailor-made molecular electronic devices. 
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1 Introduction 

Speed of growth of molecular electronics is being ac- 
celerated more and more as it has brought together 
scientists and engineers from various disciplines. 
The reason behind this attraction is inscribed into 
its smallness of size with wonderful electronic prop- 
erties. In addition, several other properties such 
as magnetic, optical, etc., have been recognized in 
different molecules, which may be utilized in arti- 
ficially tailored devices that are not possible with 
conventional materials pQ. The concept of elec- 
tron transport which emerged first in the theoretical 
work of Aviram and Ratner in 1974 [5] has opened a 
new era in the field of nanoscience. But at that time 
any type of measurement in such a small scale was 
a long-sought goal. Study at molecular scale level 
is not a simple one as we cannot avoid the effect of 
interface to the external electrodes. However, the 
progress in the theoretical works [3] was continuing, 
which bestowed inspirations to the experimentalists 
to take such task as a challenge. Now with the ad- 
vancement in nanotechnology, it is possible to inves- 
tigate several transport properties not only through 
a group of molecules [I] but also through a single 
molecule [5j. This single molecular electronics may 
play a key role in designing nanoclcctronic circuits. 
For this we have to have a thorough understanding 
of the electronic transport processes at this molecu- 
lar scale level [6j OH Q3] . Many problems are 
yet to be solved to make this field much more reli- 
able. Therefore, the electron transport in molecular 
systems is an open area and detailed investigations 
of molecular transport are still needed. 

All these works we have referred above are re- 
lated to two-terminal electron transports. We can 
also analyze various transport phenomena of a 
multi-terminal system, which was first addressed 
by Buttiker [12]. The Biittiker formalism, which is 
an extension of the Landauer two-terminal conduc- 
tance formula, is a very simple and elegant way to 
divulge the transport mechanism in terms of var- 
ious transmission probabilities. There are several 



pioneering works p~3l [14] HH US [17] based on this 
formalism, which are very interesting from the the- 
oretical as well as experimental point of view. 

Several ab initio methods [HI QH [20] [21] are 
there which may be utilized to study electron trans- 
port properties through molecular junctions. At the 
same time, tight-binding model has been extended 
to density functional theory (DFT) for transport 
calculations [22]. But in case of molecular systems, 
the investigations based on this theory (DFT) have 
some quantitative discrepancies compared to the 
experimental predictions. More over, these ab ini- 
tio theories are computationally very expensive. To 
avoid this we do model calculations by using a sim- 
ple tight-binding framework. 

In the present article we do a theoretical study 
of multi-terminal electron transport through a sin- 
gle phenalenyl molecule [23] [24] attached to semi- 
infinite one-dimensional (ID) metallic leads. We do 
exact numerical calculations based on single particle 
Green's function formalism [35] US] to evaluate con- 
ductance, reflection probability and current-voltage 
characteristics. Quite interestingly, we show that 
the positions where the leads are connected to the 
molecule as well as the presence of other leads have 
eloquent effects on these transport properties. More 
over, these characteristics are also influenced signif- 
icantly by the molecule-to-lead coupling strengths. 
These aspects can be utilized in designing nano- 
clcctronic devices. 

We organize the paper as follows. With a brief 
introduction (Section 1), in Section 2 we describe 
our model and the theoretical background. Results 
are analyzed in Section 3. Finally we conclude our 
results in Section 4. 

2 Model and a view of theo- 
retical formulation 

In this section we focus our attention on the sys- 
tems where a single phenalenyl molecule is attached 
symmetrically or asymmetrically to semi-infinite ID 



2 



metallic leads through thiol (SH bond) groups. The 
models arc shown schematically in Figs. Q] [2] and [3] 
where, the number of leads attached to the molecule 
is 2, 3 and 4, respectively. To evaluate the conduc- 
tance (g) and current (I) through this single molec- 
ular system we adopt the Green's function tech- 
nique. For this, first we define the Green's function 
for the whole system as, 

G=(E-H)- 1 (1) 

where, E = e + irj with rj arbitrarily very small 
number which can be set as zero in the limiting ap- 
proximation, e is the injecting electron energy. H is 
the Hamiltonian of the entire system which is of in- 




Figurc 1: (Color online). Two-terminal quantum 
system. A phenalenyl molecule is attached sym- 
metrically to two semi-infinite ID metallic leads, 
viz, Lead-1 and Lead-2 through thiol (SH bond) 
groups in the chemisorption technique where sulfur 
(S) atoms reside and hydrogen (H) atoms remove. 
The filled yellow circles correspond to the location 
of S atoms. 

finite dimension. So, the above equation deals with 
the inversion of an infinite dimensional matrix cor- 
responding to the system consisting of a finite size 
molecule and semi-infinite leads. However, the full 
Hamiltonian can be partitioned into sub-matrices 
that correspond to the individual sub-systems like, 

N 

h = h m + J2(h p + h pM + hI m ) (2) 
p=l 

where, Hm and H p are the Hamiltonians of the 
molecule and lead-p, respectively. N is the number 
of leads to which the molecule is attached. H p m 
represents the coupling matrix that will be non- 



zero only for the adjacent points in the molecular 
system (molecule with sulfur atoms) and the lead- 
p. Here all the leads are treated on an equal foot- 
ing. Within the non-interacting picture, the tight- 
binding Hamiltonian of the molecular system can 
be manifested as, 

i <ij> 

where, is the on-site energy, t is the nearest- 
neighbor hopping integral and cj (cj) is the cre- 
ation (annihilation) operator of an electron at the 
site i. Each lead can be described by using a sim- 




Figure 2: (Color online). Three-terminal quantum 
system. A phenalenyl molecule is attached asym- 
metrically to three semi-infinite ID metallic leads, 
namely, Lead-1, Lead-2 and Lead-3 through sulfur 
(S) atoms. 

ilar kind of tight-binding Hamiltonian, as given in 
Eq. [3l characterized by two parameters eo, the on- 
site potential and to, the nearest-neighbor hopping 
integral. 

Following the partition of the Hamiltonian, the 
Green's function can also be partitioned into sub- 
matrices and the effective Green's function for the 
molecular system can be indited (using Lowdin's 
partitioning technique [371121]) as, 

G M = - H M - J2 s ^ ( 4 ) 
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where, E p is the self-energy due to the coupling of 
the molecular system to the lead-p. It is straightfor- 
ward to obtain an explicit expression for self-energy 
corresponding to lead-p, 

s p = H pM G p H pM (5) 

where, G p = (E — Hp) -1 is the Green's function of 
lead-p. All the coupling information are inscribed 
into this self-energy expression. Once the Green's 
function is established, the coupling function T p can 
be easily obtained from the equation [29j [30] , 

r p (E) = i[z;(E)-z a p (E)} (6) 

where, the advanced self-energy S° is the Hcrmitian 
conjugate of the retarded self-energy EL Thus, we 
can write, 

r p = -2Im(£L) (7) 

In order to evaluate the conductance for the multi- 
terminal quantum system, we use the Biittiker for- 
malism }29j , valid at much low temperature and bias 
voltage, in the form, 

2e 2 

9pq = ~ T vi ( 8 ) 

where, T pq is the transmission probability of an elec- 
tron across the molecular system from the lead-p to 
lead-q and it is related to the reflection probability 
by the equation, 

R PP + E T op = 1 ( 9 ) 

which is obtained from the condition of current con- 
servation [31 j - Now, this transmission probability 
can be expressed in terms of the effective Green's 
function of the molecular system and molecule-to- 
lead coupling as, 

Tpq = Tr \rpG\[T q G a M ] (10) 

where, G r M and G^ are the retarded and advanced 
Green's functions of the molecular system, respec- 



tively. T p and T q represent the couplings of the 
molecule to the lead-p and lead-q, respectively 
Since the coupling matrix H p m is non-zero only 
for the adjacent points, n and m, the transmission 
probability becomes [32], 

T pq = 4 A p (nn) A q (mm) \ G M {nm) | 2 (11) 

where, A p (nn) = (n|A p |n), A q (mm) = (?7i|A 9 |m), 
GM(nm) = (n\GM\m) and A p , A g are the imagi- 




Figure 3: (Color online). Four-terminal quantum 
system. A phenalenyl molecule is attached asym- 
metrically to four semi-infinite ID metallic leads, 
viz, Lead-1, Lead-2, Lead-3 and Lead-4 through sul- 
fur (S) atoms. 

nary parts of E p and E g , respectively. 

In case of two-terminal system [33], Eq. [8] be- 
comes quite simpler like, 

2e 2 

.9 = -^T (12) 

and accordingly, the reflection probability becomes 
R = 1 — T. For the two-terminal quantum system 
the above expression (Eq. [T2")) is the so-called Lan- 
dauer conductance formula. 

The current I p passing through the lead-p can be 



4 



obtained from the following expression [29], 

DO 

J p = IT E / T m(E) lf P (E) - f q (E)] dE (13) 

Q 

H — oo 

where, / p ( 9 ) = f (E — /-%>(<?)) ^ s the F erm i distribu- 
tion function with the chemical potential = 
Ep ± eV p ( q )/2. Ep is the equilibrium Fermi energy. 
Throughout this calculation, we assume that the 
entire voltage is dropped across the molecule-lead 
interfaces as this assumption introduces a minimal 
effect on the behavior of the I-V characteristics. In 
this article we set c = h = e = 1 for the sake of 
simplicity. 

3 Numerical results and dis- 
cussion 

In order to illustrate the results, let us first men- 
tion the values of different parameters used in our 
numerical calculations. All the on-site energies of 
molecule, sulfur atom and leads are set to zero, 
while the nearest-neighbor hopping strengths are 
fixed at 3 for both the molecule (t) and leads 
(to). But the values of molecule-to-lead coupling 
strengths (t p for lead-p) are different from the value 
assigned for t and to. Based on the molecular cou- 
pling strength, we analyze our results in two dis- 
tinct regimes. One is the weak-coupling regime and 
the other is the strong-coupling regime. In the first 
case, t p « t and we set r p = 0.5. In the second 
case, t p ~ t and for this regime we fix t p = 2.5. 
Here we consider that r p 's are identical for all the 
leads p and set the equilibrium Fermi energy Ep at 
0. 

3.1 Conductance-energy characteris- 
tics 

In the forthcoming sub-sections we present the char- 
acteristic properties of electron transport for two-, 
three- and four-terminal molecular systems where 



the molecules are attached to leads via thiol-linking 
groups (SH bond). In experiments, leads are gener- 
ally designed from gold (Au) and the thiol groups 
are linked by using chemisorption technique [34] 
where hydrogen (H) atoms remove and sulfur (S) 
atoms survive. 

3.1.1 Two-terminal conductance 

In Fig. [4] we present the variation of two-terminal 
conductance g (red curves) and reflection probabil- 
ity R (green curves) with injecting electron energy 
E. The results in the weak molecule-to-lead cou- 
pling limit are shown in (a) and (c), while (b) and 
(d) represent the same for the strong molecule-to- 
lead coupling limit. In the weak-coupling regime, 



(a) 







(c) 




Figure 4: (Color online). Conductance g (red lines) 
and reflection probability R (green lines) as func- 
tions of energy E for two-terminal molecular sys- 
tem, (a) and (c) represent the results for the weak- 
coupling limit, while (b) and (d) correspond to the 
same for the strong-coupling limit. 

the presence of sharp resonant peaks indicates that 
electron transmission occurs at some typical en- 
ergy values, while for all other energies conductance 
vanishes (see Fig. H{a)). All these resonant peaks 
are associated with the energy eigenvalues of the 
molecular system, and therefore, we predict that 
conductance spectrum is a fingerprint of the elec- 
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tronic structure of the system. Most of the resonant 
peaks approach the value 2, the maximum value of 
g following the Landauer conductance formula (see 
Eq. [12]) and hence T goes to unity at these reso- 
nances indicating ballistic transmission through the 
molecular wire. But the behavior of conductance 
spectrum changes in case of the strong molecule-to- 
lead coupling limit. Width of each resonant peak 
becomes larger and larger as we increase gradu- 
ally the molecule-to-lead coupling strength and for 
a large molecular coupling, we have the situation 
where electron transmission takes place for the en- 
tire energy range (for illustration, see Fig. HJb)). 
The effect of such broadening comes from the imag- 
inary parts of the self-energies [29] . All these phe- 
nomena emphasize that fine tuning in the energy 
scale is necessary as long as the coupling strength 
is much weak, while it is not required in case of the 
strong molecule-to-lead coupling limit. 

This scenario is just inverted in case of reflection 
probability R. In the weak-coupling regime, sharp 
dips appear (Fig. Etc)) for some particular energy 
values where conductance shows resonant peaks, as 
R follows the simple relation R = 1 — T for the 
two-terminal molecular system. For all other en- 
ergy values, R = 1 indicating no transmission of 
electron across the molecule. The effect of molecule- 
to-lead coupling is exactly similar to that for con- 
ductance spectrum. In the strong-coupling regime 
(see Fig. Bid)), the reflection probability no longer 
reach the maximum value (1) for the entire energy 
range. 



tra. But the point is that for the three-terminal 
system most of the resonant peaks do not reach the 
value 2. From the conductance spectra it is clear 
that the heights are much reduced compared to the 
two-terminal case. This is solely due to the effect of 
quantum interference of the electronic waves pass- 
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Fi gure 5: (Color online). Conductance g P q 
(red curves) and reflection probability R pp (green 
curves) as functions of energy E for three-terminal 
molecular system. All the results are presented only 
for the weak-coupling regime. 



3.1.2 Three-terminal conductance 

To illustrate the results for the three-terminal quan- 
tum system, constructed by attaching three leads to 
the molecule (see Fig. [2]) , let us start by referring to 
Fig.[5]where the first column shows the conductance 
spectra g pq (from lead-p to lead-q) and the second 
column presents the nature of reflection probability 
R pp . Similar to the two-terminal molecular system, 
some sharp peaks appear in the conductance spec- 



ing through different arms of the molecular rings. 
In this three-terminal system, three leads are at- 
tached asymmetrically to the molecule at three dif- 
ferent locations, which provide different path ways 
for electron transmission between the leads. This 
introduces anomalous features in the conductance 
spectra as illustrated in Figs. EJa), (b) and (c). 
More over, in this three-terminal case, nature of 
variations of reflection probabilities is not so sim- 
ple as we get in the case of two-terminal molecular 
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system. Here, it is not necessary that R pp shows 
dips or peaks where g pq has peaks or dips since 
it (Rpp) depends on the combined effect of T pq s 
obeying the expression given in Eq. |9] Another im- 
portant point we like to mention here is that, one 
can easily find T qp for any two leads q and p if T pq 
is known since the relation T pq = T qp holds true 
following the time-reversal symmetry. The effect 
of molecule-to-lcad coupling strength is identical to 
the case of two-terminal conductance and therefore, 
we do not show those results further. 



andlH respectively. The results are plotted in Fig. [6] 
considering the weak molecule-to-lead coupling. In 
each figure g pq and g qp are plotted and they are 
superposed to each other due to their symmetry. 
Sharp resonant peaks of different heights for some 
particular energies appear in the conductance spec- 
tra similar to the two- or three-terminal conduc- 
tance spectra. From the conductance spectra, the 
effect of quantum interference associated with the 
molecule-to-lead interface geometry is well under- 
stood. Also the dependence of electron transport 



3.1.3 Four-terminal conductance 

The transport properties of the four-terminal 
molecular system is also described by investigat- 
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Figure 6: (Color online). Four-terminal conduc- 
tance g pq as functions of energy E in the limit of 
weak molecular coupling. 

ing various conductances g pq and reflection prob- 
abilities R pp , which are obtained from the Eqs. O 
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Figure 7: (Color online). Reflection probability R pp 
as a function of energy E for four-terminal quan- 
tum system. These results correspond to the weak- 
coupling regime. 

on molecular coupling strength is exactly similar 
to that in Fig. H] and accordingly, the results are 
not given further. Comparing the results given in 
Figs. [5] and |5] we observe that, the conductances for 
two particular leads located at the same positions of 
the molecule exhibit completely different features in 
presence of the other leads. For instance, gis of the 
four-terminal system and 1712 of the three-terminal 
system show different spectra though these two con- 
ductances are evaluated for the two leads located at 
the same positions of the molecule. Similarly, care- 
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ful observation depicts that pathways between the 
lcad-2 and lead-3 of the three-terminal system are 
exactly similar to that between lead- 2 and lead-4 of 
the four-terminal system. In spite of this the corre- 
sponding spectra i.e., 523 of Fig. [5] and #24 of Fig. [6] 
are different from each other due to the presence of 
the fourth lead. The effect of the leads are incorpo- 
rated into the self-energies which lead to the change 
in conductance spectra. Thus, we can say that the 
presence of additional leads has a pronounced ef- 
fect on the electron transport properties. In R-E 
spectra given in Fig. [7J we get various dips at dif- 
ferent energies depending on the combined effect 
of the transmission probabilities (Eq. [5]) , similar to 
the case of three-terminal system. All the reflec- 
tion probabilities are calculated only for the limit 
of weak molecular coupling. Exactly similar fea- 
ture, except the broadening, will be observed in the 
case of strong-coupling. 

3.2 Current-voltage characteristics 

The scenario of electron transport through molecu- 
lar junction becomes much more transparent when 
we discuss the current- voltage (/-V) characteristics, 
where the current is evaluated by integrating the 
transmission probability T using Eq. [T3] The na- 
ture of the variation of transmission probability is 
exactly similar to that of conductance spectra ex- 
cept the factor 2 as we have assumed e = h = 1 in 
the Landauer conductance formula (Eq. [i"2j) . Here, 
we discuss the I-V characteristics for two-, thrce- 
and four-terminal molecular systems separately in 
the following sub-sections. 

3.2.1 Two-terminal molecular system 

As an illustration, we display the I-V characteris- 
tics for the two-terminal system in Fig. 13 where 
(a) and (b) correspond to the results for the weak 
and strong molecule-to-lead coupling limits, respec- 
tively. For this two-terminal case, current can be 



expressed mathematically as follows, 
1 = 9 {Vi - V 2 ) 

= gv 12 (14) 

where, V12 is the voltage difference between the 
lead-1 and lead-2. In the case of two-terminal 
molecular system, we have attached two leads to the 
molecule and their chemical potentials are changed 
as the bias voltage is applied. With the increase of 
the voltage, the gap increases more and more and 
eventually crosses molecular energy levels one after 
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Figure 8: (Color online). Current I as a function of 
applied bias voltage V for two-terminal molecular 
system, (a) Weak-coupling limit and (b) strong- 
coupling limit. 

another. Accordingly, current channels are opened 
up and jumps in the I-V curve appear. This pro- 
vides staircase-like structure in the current-voltage 
spectrum and only for the case of weak molecular 
coupling these sharp steps appear. But this feature 
i.e., step- like changes gradually towards continuous 
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nature with the increase of molecule-to-lead cou- 
pling strength. In addition to that, current ampli- 
tude becomes much higher compared to the case 
of weak-coupling limit. By noting the area under 
the curve of Fig. 0Jb) the reason behind this en- 
hancement of current amplitude is clearly under- 
stood. Thus, it can be manifested that molecule- 
to-lead coupling strength has a significant influence 
on molecular transport. 

3.2.2 Three-terminal molecular system 

Now, we describe the current-voltage characteris- 
tics for the three-terminal system where we find the 
current I p in lead-p by integrating the transmission 
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Figure 9: (Color online). I\ as a function of V13 
(= Vi — V3) for the three-terminal molecular system 
in the limit of strong-coupling, keeping V12 (= Vi — 
V2) as a constant. The orange, green and magenta 
colors correspond to V12 = 1,2 and 3, respectively. 

function T pq considering the effects of other termi- 
nals. To be more precise, we write the current ex- 
pressions for three different leads as follows, 

h = 912 {Vi - V 2 ) + gi 3 (Vi - V 3 ) 

= 312V12 + 513V13 (15) 
h = 921 { V 2 - Vi ) + .g 23 (V 2 - V 3 ) 

= 521^21+023^23 (16) 
^3 = 531 (V3 -Vi)+ 932 (V 3 - V 2 ) 

= 531 V31 + 532 V32 (17) 



where, V pq — (V p — V q ) is the voltage difference be- 
tween the two leads named as lead-p and lead-q. 

In Fig. |9j we plot I\ for the lead-1 as a function 
of V\ 3 for different fixed values of V12 in the strong 
molecule-to-lead coupling limit. The orange, green 
and magenta curves represent the currents for V12 = 
1, 2 and 3, respectively. It is clear from the figure 
that for a constant value of V12, the moment we 
switch on the bias voltage between lead-1 and lead- 
3, current rises to a large value. Then, for a wide 
range of V13, it (Ii) slowly increases with the rise of 
V13 and finally the rate of increment of the current 




V21 



Figure 10: (Color online). I2 as a function of V21 
(= V 2 — Vi) for the three-terminal molecular system 
in the strong-coupling regime, keeping V23 (= V2 — 
V3) as a constant. The blue, pink and green lines 
correspond to V23 = 1,2 and 3, respectively. 

gets enhanced when V13 is quite high. This behavior 
i.e., the rate of increment of the current with bias 
voltage solely depends on the positions of resonant 
peaks in the g-E spectrum. For a particular value of 
V\ 3 , 1 1 increases as we increase V12 which is clearly 
visible from three different curves plotted in Fig. [SJ 

The variation of current I2 through lead-2 as a 
function of V21 is shown in Fig. 1101 keeping the 
voltage V23 as a constant. The currents are eval- 
uated in the strong-coupling limit, where the blue, 
pink and green lines correspond to V23 = 1, 2 and 
3, respectively. 

In a similar fashion in Fig. [TT] we display I 3 - 
V31 characteristics for different fixed values of V32 
considering the case of strong-coupling limit. The 



9 



magenta, green and orange curves correspond to 
V32 = 1, 3 and 5, respectively. 

The characteristic features of the currents Ii, 
I 2 and ^3 passing through three different leads 
are quite analogous to each other. Depending on 
the conductance-energy spectra, we get different 
current amplitudes for three different leads which 
are clearly observed from the results presented in 
Figs. [9l[T0l and [TTl All these currents are computed 
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Figure 11: (Color online). ^3 as a function of V31 (= 
V3 — Vi) for the three-terminal molecular system in 
the strong molecule-to-lead coupling limit, keeping 
V32 as a constant. The magenta, green and orange 
colors correspond to V32 = 1,3 and 5, respectively. 

only for the strong-coupling limit. We can also de- 
termine the currents for the limit of weak-coupling 
and in that case we will get sharp step-like features 
as a function of bias voltage with much reduced am- 
plitude compared to the strong-coupling case. The 
origin of step- like behavior in current is clearly men- 
tioned in the case of two-terminal molecular system 
(Sec. 3.2.1). 

From the above current expressions (Eqs. [T5l [TO] 
and I17|) we see that the current in anyone lead is 
related to two potential functions. For instance in 
Eq. 1151 there are two parameters like V12 and V13 . 
Keeping V\ 2 as a constant we plot the current 1\ in 
terms of V13 (see Fig. [9]). At the same time, we can 
also draw the h-V\ 2 characteristics, considering V13 
as a constant. Both for these two cases, the char- 
acteristic features are quite similar. This argument 
is also valid for the other two current expressions 



(Eqs. H6]and [H]). 

All the above features of current-voltage char- 
acteristics in this three-terminal molecular system 
clearly support the basic features of a traditional 
macroscopic transistor. Thus we can predict that 
the three-terminal molecular system may be utilized 
to design a nano-scale molecular transistor. 

3.2.3 Four-terminal molecular system 

Finally, we focus our attention on the current- 
voltage characteristics for the four-terminal molcc- 
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Figure 12: (Color online). 1\ as a function of V13 for 
the four-terminal molecular system in the strong- 
coupling limit, keeping V\ 2 and V14 as constant. 
The magenta, green and orange colors correspond 
to V\2 = V14 = 1,3 and 5, respectively. 

ular system. In this molecular system, the current 
expressions for the four different leads arc as fol- 
lows, 

h = 9i2 (Vi - V 2 ) + 913 (Vi - V 3 ) + gu (Vi - V 4 ) 

= gnVi2 + 013 Vl3 + 014 V14 (18) 

h = 921 (V 2 -Vi) + 323 {V 2 -Vk)+ S24 {V a - V 4 ) 

= 92 1 V 2 1 + 323^23 + 324 V24 (19) 

h = 931 (Va - V x ) + 932 {V 3 - V 2 ) + g 3 4 {V a - V 4 ) 

= 331^31 + 332 V32 + 334 V34 (20) 

h = gii{Vi-Vx) + g i2 {Vi-V 2 ) + g i a{V i -Va) 
= gaVii + 342^42 + 343 V43 (21) 

where, V pq is the voltage difference between the 
lcad-p and lead-q. 
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As representative examples, in Fig. [12] we plot the 
current in lead-1 as a function of V13 keeping 
V12 and Vu as constant. The results are computed 
for the strong-coupling limit, where the magenta, 
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Figure 13: (Color online). I2 as a function of V24 
(= V2 — V4) for the four-terminal molecular system 
in the limit of strong-coupling, keeping V21 and V23 
as constant. The green, pink and dark-blue lines 
correspond to V21 = V23 = 2, 4 and 6, respectively. 

green and orange curves correspond to V12 = V14 = 
1, 3 and 5, respectively. The variation of current 
1\ as a function of V13 in this four-terminal molec- 
ular system is quite similar to that as presented in 
the case of three-terminal system (Fig. For a 




Figure 14: (Color online). Is as a function of V31 
for the four-terminal molecular system for the case 
of strong-coupling limit, considering V32 and V34 as 
constant. The green, magenta and dark-blue curves 
correspond to V32 = V34 = 1, 3 and 5, respectively. 

particular value of V13, here also the current am- 
plitude gets increased with the rise of V12 and V14. 
But quite significantly we observe that, for a par- 



ticular value of V13 , the current I\ passing through 
the lead-1 of four-terminal system acquires much 
higher amplitude compared to the three-terminal 
system (see Figs. [9] and Q2]). The reason behind 
this enhancement of current amplitude is explained 
as follows. From Eq. [18] we see that I\ contains 
three additive terms where the contributions come 
from other leads, while in Eq. [15] there are two ad- 
ditive terms. The additional term appears in Eq. 1181 
is due to the presence of fourth terminal which is 
responsible for the larger current in four-terminal 
system compared to the three-terminal one. 

In case of the current I2 through lead-2 we show 
the variation with respect to V24, keeping V21 and 
V23 fixed to a particular value as presented in 




Figure 15: (Color online). I4 as a function of V43 
(= V4 — V3) for the four-terminal molecular system, 
keeping V41 and V42 as constant. The currents are 
evaluated in the strong-coupling limit, where the 
green, dark-blue and pink lines correspond to V41 = 
V42 = 1, 5 and 9, respectively 

Fig. 1131 The currents are determined in the strong- 
coupling regime, where the green, pink and dark- 
blue lines correspond to V21 = V23 = 2, 4 and 6, 
respectively. 

Similarly, in Fig.Q3]we display /3-V31 characteris- 
tics considering V32 and V34 as constants in the limit 
of strong-coupling. The green, magenta and dark- 
blue lines represent the currents for V32 = V34 = 1, 
3 and 5, respectively. 

At the end, in Fig. [15] we show the variation of 
current I 4 passing through lead-4 as a function of 
V43 in the limit of strong molecular coupling, when 



11 



Vn and V42 are kept as constants. The green, dark- 
blue and pink curves correspond to the currents for 
V41 = V42 = 1, 5 and 9, respectively. 

Similar to the case of three-terminal molecular 
system, for this four-terminal case, we can also plot 
the current through any lead-p in aspect of the 
other potential functions as given in Eqs. [THl HSl 
l20l and [2TJ For all these cases, the characteristic 
features are very much similar those are presented 
in Figs. H21 HH and [15l and hence, we do not 
re-plot the results further. 

4 Closing remarks 

In the present paper we have used a parametric ap- 
proach to study multi-terminal electron transport 
through a single phenalenyl molecule. Using a sim- 
ple tight-binding framework we have performed all 
the numerical calculations through single particle 
Green's function formalism. The basic features of 
electron transport in this molecular system are ex- 
plored by investigating the multi-terminal conduc- 
tance, reflection probability and current. Following 
a detailed description of electron transport in two- 
terminal quantum system, we have revealed the es- 
sential features of electron transport in the three- 
and four-terminal quantum systems separately. 

Our clear investigation predicts that the electron 
transport in multi-terminal molecular system signif- 
icantly depends on (a) the molecule-lead interface 
geometry, (b) the presence of other leads and (c) the 
strength of molecular coupling to the side attached 
leads. The unique characteristics of this phenalenyl 
molecule with a very small size has enhanced the 
importance of the present article. Our parametric 
study provides several significant features to reveal 
electron transport through any complicated multi- 
terminal quantum system. 

In the present paper we have done all the cal- 
culations by ignoring the effects of temperature, 
electron-electron correlation, etc. We need further 
study by incorporating all these effects. 
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